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Abstract. We study generating functions in the context of Rota-Baxter algebras. We 
show that exponential generating functions can be naturally viewed in a very special case 
of complete free commutative Rota-Baxter algebras. This allows us to use free Rota-Baxter 
algebras to give a broad class of algebraic structures in which generalizations of generating 
functions can be studied. We generalize the product formula and composition formula 
for exponential power series. We also give generating functions both for known number 
families such as Stirling numbers of the second kind and partition numbers, and for new 
number families such as those from not necessarily disjoint partitions and partitions of 
multisets. 
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The importance of generating functions is well-known [30, 31]. We view generating 
functions in the framework of Rota-Baxter algebra whose importance in combinatorics was 
already noted by Rota |27], ^9|. This viewpoint leads to a broad generalization of exponential 



generating functions. 

There are many types of generating functions, such as the power series (or ordinary) gen- 
erating functions f(z) = a n z n , exponential generating functions f(z) = Yl a n^-j and 

n>0 n>0 

multi-variable generating functions. For some number sequences (a n )„>o, their exponential 
generating functions have good expressions that are not available if power series generating 
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functions are used. Such examples include the Bernoulli numbers and Bell numbers. Fur- 
ther, multi-variable generating functions encode multi-indexed number families that cannot 
be represented by one variable generating functions. 

One way to view the power series generating functions and exponential generating func- 
tions in the same framework is that they each give a way to encode a number sequence 
a n ,n > as the coefficients of a linear combination with respect to a basis of the power 
series algebra The basis is z n ,n > for power series generating functions and is 

the divided powers z n /n\,n > for exponential generating functions. The distinction is 
the two different bases for the same power series ring. But it is also beneficial to view the 
difference externally: V). = z k , k > is the standard basis of the algebra 

A := Y\_^ v k, with componentwise product v m v n = v m+n , 

k>0 

while Wk = z k /k\ is the standard basis of the divided power algebra 



H> := J^Mu^., with componentwise product 

fc>0 




We can take this point of view further and consider the following general framework for 
generating functions: A complete filtered M-algebra is a M-algebra A with ideals A n , n > 0, 
such that A m A n G A m j rn and A is complete with respect to the metric on A induced by 
A n . In other words, the natural map 

A -)> ]imA/A n 

is bijective. Let 11 := {uj,j G J}, be a basis of A that is compatible with its filtration 
in the sense that U R Aj. is a basis of A/~, k > 0. Then a U-generating function of 
a family of numbers dj G K,j G J, is the element X]jeJ a i M i * n ^" ^ ^ ms con t ex t, a 
power series generating function is a U-generating function when XL is taken to be the basis 
Vk = {z k ,k > 0} in the complete filtered algebra 

R[[z\] = Y[Rz k = Y[Rv k = A 

k>0 k>0 

and an exponential generating function is a U-generating function when U is taken to be 
the basis {w k = z k /k\, k > 0} in the complete filtered algebra 

k 

r[[z]) = n M f[ - n Mujfc=:B - 

fc>0 ' fc>0 

In both cases, the complete filtration on R[[z}} is given by the ideals z n R[[2:]] which is also 
UkA (resp. Vk'B). 

Of course such a formal definition in such generality is of little use unless 

(1) it can be naturally related to the ordinary generating functions or exponential gen- 
erating functions; 

(2) it is useful in the study of number sequences and number families. 

We will show that free Rota-Baxter algebras do give a generalization that satisfies these 
conditions. In Section [2| we review the construction of free commutative Rota-Baxter al- 
gebras and show that their completions give a large class of complete filtered algebras. In 
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Section [| we show that the simplest case of such complete free Rota-Baxter algebras gives 
the exponential generating functions. We then focus on a "twisted" variation of this case, 
namely on the A-exponential generating functions. We generalize the classical Product For- 
mula and Composition Formula for exponential generating functions to the twisted case. 
As applications, we give generating functions for Stirling numbers and for enumerations 
of not necessarily disjoint partitions. We consider another instance of complete free Rota- 
Baxter algebras in Section [|. By using g-series and applying the stufHe technique from the 
study of multiple zeta values in number theory, we show that elements of these algebras 
provide interesting generating functions for both number sequences such as factorials of 
square numbers, triangular numbers, pentagonal numbers, and multi-indexed families of 
numbers such as those from partitions of multisets. 



2. Free commutative complete Rota-Baxter algebras 

In this section, we provide the algebraic framework for our extension of exponential 
generating functions. Let A 6 R be a constant. A Rota— Baxter algebra of weight A is a 
pair (R, P) where R is a unitary k-algebra and P : R — > R is a linear operator such that 

(1) P(x)P(y) = P(xP(y)) + P(P(x)y) + XP(xy), 

for any x, y G R. Often 9 = —A is used, especially in the physics literature. Let (R, P) and 
(R', P') be Rota-Baxter algebras of weight A. A Rota-Baxter algebra homomorphism 
/ from (R, P) to (R', P') is an algebra homomorphism / : R — > R' such that / o P = P' o f. 
The study of Rota-Baxter algebras was started by G. Baxter [§] in 1960 and was popu- 



larized largely by Rota and his school in the 1960s and 70s |27|, g8] and again in 1990s [^9 
In the recently years, there have been several interesting developments of Rota-Baxter 
algebras in theoretical physics and mathematics, including quantum field theory, Yang- 
Baxter equations, shuffle products, operads, Hopf algebras, combinatorics and number the- 
ory [@, g 0, |TJ, |2TJ . See jT7| for further details. 



2.1. Free commutative Rota-Baxter algebras. We recall the construction of free com- 



mutative Rota-Baxter algebras in terms of mixable shuffles. See []18 , 19[ [L6[ for details. 

Given a commutative algebra k which will often be taken as IR, A G k, and a k-algebra 
A, the free commutative Rota-Baxter k-algebra on A is defined to be a Rota-Baxter k- 
algebra ( UI^xi^A) , Pa) together with a k-algebra homomorphism ja '■ A — > UI^^A) with 
the property that, for any Rota-Baxter k-algebra (R, P) and any k-algebra homomorphism 
/ : A — >■ R, there is a unique Rota-Baxter k-algebra homomorphism / : ( UI^^A) , Pa) — > 
(R, P) such that ja ° / = / as k-algebra homomorphisms. 

One realization of this free commutative Rota-Baxter algebra is given by the mixable shuf- 
fle Rota-Baxter algebra. The mixable shuffle Rota-Baxter algebra is a pair ( IIIk^A), Pa), 
where UI^xiA) is a k-algebra in which 

• the k-module structure is given by the direct sum 

oo 

A® n , where A® n = A ® k . . . ® k A; 

11 n— factors 
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• the multiplication is given by the augmented mixable shuffle product o, recursively 
defined on A® m g A® n by 

ao o (bo g h g> . . . <g> 6 n ) = a &o ® &i <8> • • • g &n, 
(a g ai g . . . g a m ) o 6 = a-o&o g ai g . . . g a m , a», G A® 1 = A, 
and 

(2) (a g ai g • • • g a m ) o (6 g &i <E> • • • g 6 n ) 

= (a &o) g ((ai g • • • g a m ) o (1 <g 6i g . . . <g> 6 n )) 
+ (a 6 ) g ((1 g ai g . . . g a m ) O (b x g . . . g b n )) 
+\a b <g ((oi <g) . . . (g) a m ) o (&i (g) . . . (g) fe n )) , a i} bj G A 
with the convention that 

a O (1 <g b) = a <g b, (1 g) a) o b = b g) a, aob = ab, for a, 5 G A 
The Baxter operator Pa is defined by 

Pa(oi g) . . . (g> a m ) = 1 g) ai g . . . !g> a m , ai g . . . g a m G A® m , m > 1. 

Since the mixable shuffle product is compatible with the product on A, we will often suppress 
the notation o. We will also suppress k and A from Ul^x(A) when there is not danger of 
confusion. 

Note that assuming Pa is a Rota-Baxter operator and thus satisfies Eq. ([I]), then Eq. (0) 
follows. For example, we have 

(a g a x g a 2 )(& g fej) 
= a &o g ((ai g a 2 )(l g &i) + &i(l g a x g a 2 ) + A(a x g a 2 )&i) 
= ao^o g (a-i g (0.2(1 g 61) + (&i(l g 0-2)) + Aa 2 6i) + 61 g a x g a 2 + Aai&i g a 2 ) 
= a &o g (ai g a 2 g 61 + a\ g 61 g a 2 + Aai g a 2 6i + 61 g a x g a 2 + \a\b\ g a 2 ) . 

Theorem 2.1. ([[l|, Theorem 4.1]) For any k-algebra A, ( Ul(A),P A ), together with the 
natural embedding ]a '■ A — )■ IIl(A) ; zs a /ree Baxter h-algebra on A (of weight X) in 
the sense that the triple ( Ul(A) , P A , j a) satisfies the following universal property: For any 
Baxter k-algebra (R, P) and any It-algebra map (p : A — >■ R, there exists a unique Baxter 
k-algebra homomorphism (p : ( HI (A), Pa) — > (R,P) such that the diagram 



m(A) 




commutes. 

Alternatively, Hl(A) can defined to be the tensor product algebra A® IH + (A) where the 
multiplication on IH + (A) = © fc>1 A® k is given in the explicit form by the mixable shuffle 
product |0|] and in the recursive form by a generalization |2(| of the quasi-shuffie 



algebra defined by Hoffman | in the study of multiple zeta values. It is also shown that 



HI + (A) is the free commutative tridendriform algebra |26[ 
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Quasi- shuffle product is also known as harmonic product |22| and coincides with the 



stufffe product || |f in the study of multiple zeta values. Variations of the stuffle product 
have also appeared in [[7], [TTj] . It is shown || to be the same as the mixable shuffle prod- 



uct [Tj| which is also called overlapping shuffles |JT| and generalized shuffles [|TJ[], and 



can be interpreted in terms of Delannoy paths |], |12|, |26 



2.2. Complete free Rota-Baxter algebras. For a given commutative algebra A, it is 
easy to see that the submodules 



Fil fc Hl(A) := A®\ i > 



i>k 



of III (A) are ideals of HI (A). They are in fact Rota-Baxter ideals of HI (A) in the sense 
that P A (Fi\ k IH(A)) C Fil fc m(A). Further, n fc > Fil fc Hi (A) = 0. Thus 

1h(A) := lim Hl(A)/Fil fc Hl(A) JJ A® fc 

fc>0 

is a complete filtered algebra and contains HI (A) as a subalgebra. It coincides with the 
complete free commutative Rota-Baxter algebra defined in 

3. Generating functions from Rota-Baxter algebras 

We first interpret exponential generating functions in terms of the special case of free 
commutative Rota-Baxter algebra when the ring A is the base field R and when the weight 
A G R is zero. We then consider the nonzero weight case (still taking A = R) and generalize 
to this case the Product Formula and Composition Formula for exponential generating 
functions. Applications of these formulas are provide. Other instances of free Rota-Baxter 
algebras will be investigated in the next section. 

3.1. Connection with exponential power series. Let A = R. Then 

(3) HI A (A) = HI A (R) = © fc >iR 0fe = © fe >iR l k , 

where = 1 © — ■ © 1 . The augmented shuffle product in this special case is 

(fc+l)— terms 

min(m,n) 

L ( m+n—k ^ I m 

k=0 



( 4) i m oi„= »'r' Mi 



'-m+n— k- 



When A = 0, we have 



i i _ | m+n i i 



giving the divided powers. We then have HIa(R) = n/c>o-^l fc - ^his is also the cofree 
differential algebra [^4|] with the differential operator d(x n ) = x n -i,d(l) = 0. 
Denote x n = x n /n\ (divided powers). Then as an algebra, 

= ©„> R^n 
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with multiplication given by x m x n = I m+n ) x m+n . This extends to an isomorphism 

m 



R[[x]] -> m(R) = JJrI*, x k h-> l fc , fc > 0. 

fe>0 

Through this isomorphism the theory of exponential generating function is translated to a 
theory of generating functions in Hl(R), and can be generalized to JR\{A) for other algebras 
A and other weight A. We will next set up the foundation of our approach by proving the 
product formula and composition formula for A-exponential generating functions. 

3.2. Product formula for A-exponential generating functions. We first consider free 
Rota-Baxter algebra ni(R) on R of weight 1. It is given by Eq. (fj) with product given by 
Eq. (D with A = 1. 

We quote from |H| the following simple yet fundamental property of exponential gen- 
erating functions which underlies the prominent role played by these generating functions. 



See |30[ for details. For any function / : N — > R, let 

%(*) = £/( 



n) — 
n\ 

k>0 



be the exponential generating function of /. 



Proposition 3.1. (Stanley |5(J, Proposition 5.1.1]) Let j^Y be the cardinality of a finite set 
Y . Given functions f, g : N — > R, define a new function h : N — > R by the rule 

(5) K4X) = f(#s)g(#T), 

(S,T) 

where X is a finite set, and where {S,T) ranges over all weak ordered partitions of X into 
two blocks, i.e., S D T = and S U T = X . Then 

(6) E h (x) = E f (x)E g (x). 



Definition 3.2. For / : N — )■ R, we call 

(7) E XJ :=J2f(k)lk£ niA(R) 

the A-exponential generating function (or A-EGF in short) of /. When there is no 
danger of confusion, we will suppress A from the notation. 

When A = 0, this recovers the exponential generating function. We prove the following 
generalization of Proposition |3.1| . 

Proposition 3.3. Given functions f, g : N — >■ R, define a new function h : N — > R by the 

rule 

(8) M#*)= E \* {snT) f(#s)gm), 

(S,T)e7(X) 2 ,SuT=X 
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where X is a finite set and 7(X) is the power set of X . Thus in the sum (S,T) ranges 
over all ordered pairs of subsets of X such that S U T = X (but not necessarily S PI T = 0). 
Then 

(9) Ex,h = E X jE x , g g m A (R). 



Proof. We have 

E XJ E X>9 = (^/(m)l m )(^ 5 (n)l n ) 



m>0 n>0 



rri\n 

m,n>0 

min(m,n) , \ / \ 

e /(».)»(») e ^(tm:) w. 



m,n>0 i=0 

by Eq. (H). Note that I m+ ™ ' j I m j = I m+n 2 1 . So setting u = m + n — i, u\ = i,U2 

m I \ i I \ i,m—i,n—i 



m — i, -u 3 = n — i, we have u = U\ + U2 + u% and 

E XJ E x , g = J2( V A M1 ( u )/(«i + « 2 )y(tti + « 3 ))l„. 
u =0 (ui,u 2 ,M 3 )eN 3 

Now the theorem follows since, for given (ui, v,2, M3) G N 3 , I ] is the number of ways 

y u\,u2,u3 J 

of partitioning a size u set into three subsets of size «i, «2 and u%. It is also the number of 
ways of taking subsets S of size Ui + u 2 and T of size Ui + u 3 of a set X of size tti + w 2 + M3 
such that S U T = X and #(£ n T) = u v □ 



From the proof of Proposition |3]3| we have the following explicit formula for h(#X) in 
Eq. (§. 

Corollary 3.4. The numbers in Eq. (fjj is given by 



(10) h{u) 



J2 XU1 ( U )f(ui + U2)g(ui + u 3 ),u>0. 



(«i,«a,«3) G N 3 

"1 + U 2 + M3 = U 



In analogy to [3C], we have the following combinatorial significance of Proposition 3.3 



Suppose we have two types of structures, say a and j3, which can be put on a finite set 
X. We assume that the allowed structures depend only on the cardinality of X. A new 
"combined" type of structure, denoted aU x (3, can be put on X by placing structures of type 
a and /3 on subsets S and T, respectively, of X such that S U T = X (but not necessarily 
SnT — 0). If /(&) and g(fc) denote the sums of allowed structures on a k-set of type a and 
(3, respectively, then the right-hand side of Eq. (H) counts the sum of the allowed structures 
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of type aU\ /3 with a "correction factor" \#( SnJ ") that measures the overlap of S and T. In 
particular, when A = 0, then we recover the classical case in [3D . 



Example 3.5. As an application, suppose a group X of n children participates in an event 
where they can play two games a and (3. Of course a child can play either one or both 
or none of the games. Let the group of children who play game a and (3 be S and T, 
respectively. Suppose game a is a simple hit-or- missing game (such as shorting a ball). So 
there are 2* s possible outcomes. Suppose game is a competition that results a linear 
order of the participants. So there are (#T)! possible outcomes. Let h(n) = h(#X) be 
the possible outcomes of the event. Let E\j = J2 k>0 f(k)lk and E\ j9 = ^2 k>0 g(k)lk be 
the A-EGF for the outcomes of game a and (3, respectively. Then the A-EGF for the whole 
event is 

E\,h — E\jE\ >g . 

Using an inductive argument, we obtain the following generalization of Proposition 5.1.3 
in [001. 



Theorem 3.6. (Generalized Product Formula) Fix k G P with k > 2 and functions 
fi, • ■ ■ , fk : N — >• JR.. Define a new function h : N — > K. by 

(11) KW)= £ A# Tl+ - + # T -#^ u - u ^/ 1 (#T 1 )---/ fc (#T fc ), 

(Ti,-,T k ) 



or 



(12) *i(#T)= £ A^'-' Tfe )/i(#Ti)---/ fc (#T fc ), 

(Ti,-,T fc ) 

where (Ti, • ■ • , T^) ranges over all weakly ordered (not necessarily disjoint) subsets Ti, • - • , T k 
of S such that T\ U • • ■ U T k = T and 

v (T u ...,T k )= (-i) #/ #(D T 0- 

i c [k] iei 

Then 

k 

(13) E x , h = HE XJi 

i=i 

in Jn. x (R). 

Proof. We first use induction on k > 2 to prove that the fuction h : N — > M. defined by 
Eq. (0 satisfies Eq. flTj|). Since #(T X n T 2 ) = #T X + #T 2 - #(T X U T 2 ), the case of fc = 2 
is proved in Proposition |3.3| . Assume that the claim has been proved for k = n > 2 and 
consider the case of k — n + 1. Let /i, • • • , / n +i : N — > R be given. Let /i n : N — > R be 
defined such that 

E\j ln = E x j x ■ ■ -E\j n . 
Then by the induction hypothesis, 

(14) /m(#T)= Yl A^-^-^-^/U^)---/^^). 

(Ti,-,T„) 
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Since E\j = E\j ln E\j n+1 by definition, by Proposition |373| (that is, the case when k = 2) 
we have 



(15) 



/(#*) 



£ A^ + #^-#^ u ^)/ ln (#X 1 )/„ +1 (#X 2 ). 



(Xi,x 2 ) e y(x) 2 
Xi u x 2 = x 



Combining Eq. (]14|) and Eq. ([15]) we have 
/(#*) 

/ 

= ^2 \* Xi+#X2 ~* {XiuX2) 

(Xi,x 2 ) e 7{x) 2 



\ 



Xi U X 2 = X 



(Xi.i,... ,Xi, n ) eo>(Xir 

\ Xi,i U---UXi,„ =Xi 



/n.+ l(#X 2 ) 



/ 



£ A #x 1+ #x 2 - # (x 1 ux 2 ) +E?=1# x 1 , i - # (u ?=1 x li0/l(#Xlil) . . . f n{#Xl , n )f n+1 (#X 2 ) 



(Xi,i,--- ,Xi, n ,X 2 ) G ^(X)^ 1 
Xi i U ■ • • U Xi „ U X n +i = X 



since any (X ljl; ■ • • , X lin , X 2 ) G T(X) n+ such that X lj:L U- ■ -UX lin UX n+1 = X corresponds 
uniquely to a (Xi,X 2 ) G IP(X) 2 with X x U X 2 = X together with a (X^i, • • • ,X^ n ) G 
7{Xi) n with Xi,! U • • • U Xi <n = X\. Thus to complete the induction we only need to show 

n 

#X! + #x 2 - #(x x ux 2 ) + X: Pw - #(ur =1 X! 



=#x 1;1 + ■ • • #x ljn + #X 2 - #(X Xil U ■ • ■ U X 1>n U X 2 ) 

which is clear since Xi = X ltl U ■ ■ • U X l n . 

Now to finish proving the theorem, we just need to show that Eq. ([11]) agrees with 
Eq. (ID, that is, 

#T X + • • • + #T fe - #(T X U • • • U T fc ) = ^(Tx, • • • ,T k ). 
This follows from the well-known Euler characteristic for the subsets T 1; ■ ■ • ,T k : 

A: 

#U T * = E (-l) # ^ 1 #(n ie /T l ) = #T 1 + --- + #T fc -r ? (T 1 ,--- ,T k ). 

$^IC[k] 

□ 



For the differential operator g^, we have ^(ln) = ln-i, ^a(I) = 0, and the following 
result. 

n>0 n>0 

As the applications of Theorem |3.6| , we have the following results. 



Corollary 3.7. Let S be a finite set. Given functions f, g : N — > M ; define new functions 
hi, h 2 , ti3, and h^ on N as follows: 

hi(#S) = f(#S) + g{#S) 

h*(#S) = (#W(#T) + A(#£)/(#£), where #T = #S-1 
hs{#S) = /(#T), where #T = #S + 1 
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= (#5)/(#5) + A(#5)/(#r), tflfcere#T = #£ + l. 

T/ien 

-^A,?i2 = IiO-E'a,/ 
E\hz = d A (E X j) 
Exm = liod A (E x ,f). 

The proof is easy and will be omitted. 

Corollary 3.8. ( fllql ) For X G R, we have the generating function 

j oo n oo oo 



1 - (1 <g> l)u 

v 7 n=0 fc=0 n=0 fc=0 

where S(n, k), n, k > 0, are i/ie Stirling numbers of the second kind. 

Proof. Since = XT =0 (1 ® 1 ) n «"> we set ^ = E^UM*) 1 * = ( l ® and 

^a,/ = Em=o = 1 ® 1- Hence, we have 

E\ h — E ™ f = J5\ / o -Ea,/ o • • • o E'a,/ . 

v M / 

n 

For i?A,/ = 1 <8> 1, we have 

r/ s II, m = 1 



0, m^l. 
According to Theorem |3l|, we get 

(16) h(#T)= Yl A^ + - + #^#^ u - uT ")/(#T 1 )---/(#T n ), 

(Ti,-,T„) 

where (T 1; ■ • • , T n ) ranges over all weak ordered (not necessarily disjoint) subsets T 1; ■ • • ,T n 
of T such that T x U • • • U T n = T. 

Notice that the nonzero terms in the right hand side of (|i6|) are those which satisfy the 
condition = . . . = #T n = 1. Therefore, we can reformulate fll6D as follows. 

h(k)= 

Si,...,Sk 

where (S±, ■ • • , S^) ranges over all the ordered partitions of [n] such that Si fl Sj = for 
1 < i, j < and Si U • • • U Sk = [n]. That is to say, 

h(k) = k\S{n,k)X n - k . 

Hence, we get 

n 

(17) (1 ®1)» = £ fc!S(n, A;)A"^l fe . 

fc=0 

Then the corollary follows. □ 
From the above proof, we have the following generation function. 
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Corollary 3.9. For A 6 1, we have the generating function 



e (l®l)u _ 

n=0 k=0 n=0 k=0 

Setting 

oo 

E\j := In, 



n=0 



we consider the expression for E^j. 
Corollary 3.10. For A = 1, we have 



Eii = X> fc - l ) ni «- 



n=0 



Proof. By Theorem 3.6, we have 



(is) ^ = (E 1 »)* = E E ^ 

n=0 n=0 (Ti,- ,T fc ) 

where (Ti, • • ■ , Tj.) ranges over all weakly ordered (not necessarily disjoint) subsets T\, • ■ ■ 
of [n] such that T\ U • • • U T*. = [n]. Let 

/ i: [n]^{0,l}, /,(j) = ( n ' ^55' 1<^<^, 



L 0, r 4 , 

be the characteristic function of Tj, 1 < z < fc. Then 

:= {* G [A;] I fi{j) ^ 0}, 1 < j < n, 
are nonempty subsets of [k]. Conversely, given nonempty subsets Sj,l < j < n, define 

9j : [k] {0, 1}, #(*) = { ^ l<i<n. 

Then 

^:={jeM|^W^o}, i<*<fc, 

form a weakly ordered (not necessarily disjoint) subsets of [n] such that 7\ U • • • U = [n]. 
Thus the set of such weakly ordered subsets {T 1; ■ • • ,T k } is in bijection with the set of 
nonempty subsets {Si, ■ ■ ■ , S n } of [k]. Therefore in Eq. fllTf), we have 

E i = (2 fc -ir, 

(Ti,-,T fc ) 

hence the corollary. □ 

Now we give a notation that will be used in the next corollary. 

Definition 3.11. For k, £ > 1 and I < n < k£, let B(n, k, £) denote the number of /c-tuples 
(Ti, ■ ■ ■ , Tfc) of (not necessarily disjoint) subsets Ti, • ■ ■ ,T& of [n] such that TiU- • -UT^ = [n] 
and #Tj = £ for 1 < i < k. For other values of n, fc, we set B(n, k, £) = 0. 
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Corollary 3.12. For £ > 1, we have the generating function 

oo hi 



yzy £ = 5>) fc = E E xM ' nB ^ k > Q 1 *- 

k=l k=l n=l 

Proof. Set 

oo 

E X j :=le = 22 f(m)l m , where f(m) = St >m . 

m=0 

Then by Theorem |3.6| , we have 

oo 

(^) fc = E E ^ Ti+ --- + * Tk -* (TiU --- UTk) f(m)---m^^ 

n=0 (Ti,-,T k ) 
ik 

(19) =J2^ n B(n,k,£)l n , 

n=l 

giving the corollary. □ 
In Eq. (|T9"D, when £ = 1, we get Eq. ([PTD; when A = 0, since B(k£, k, £) = we obtain 

(1 ® I 91 )* = ^J(1 »!•")■ 

So we get the following result. 

£ fc=o fc=o y ' 

3.3. Composition formula for A-exponential generating functions. We now give a 



generalization of the Composition Formula [30 



Theorem 3.13. (Composition Formula) |30| Theorem 5.1.4] Given functions f : P — > 
and g : N — > R roi/i g(0) = 1, define a new function h : N — > R by 

(21) h(#S) = J2 f(#Bi)'~f(#Bk)g{k),#S>0, 

TT={B u -,B k }&l(S) 

(22) h(0) = 1, 

where the sum ranges over all partitions ir = {B±, ■ • • , Bk} of the finite set S. Then 

(23) E h (x) = E g (E f (x)). 



Now we give some definitions before we prove the main theorem for the composition rule 
in IHa(R). 

Definition 3.14. A collection {Bi, • • • , B^} of [n] is called a generalized partition (with 
distinct max) of [n] if 

(1) Bi ^ 0, 1 < i < k; 

(2) B 1 UB 2 U---UB k = [n]; 

(3) maxi?! < max5 2 < ■ ■ ■ < maxi^. 
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Let n'([n]) denote the set of generalized partitions of [n]. Let S(n,k) denote the number 
of the generalized partitions of [n] with k blocks. By convention, we put 

5(0,0) = l,S{n,0) = for n > 1, and ~S(n, k) = for k > n > 1. 

Let B(n) denote the number of the generalized partitions of [n]. That means B(n) = 
T^=iS(n,k). 

For example, S(3, 2) = 8. All the generalized partitions of [3] with 2 blocks are listed 
below. 

{{1},{2,3}}, {{1},{1,2,3}}, {{1,2},{3}}, {{1,2}, {2, 3}}, 
{{2},{1,3}}, {{2},{1,2,3}}, {{1,2}, {1,3}}, {{1, 2}, {1, 2, 3}}. 

But {{1, 3}, {2, 3}} and {{1, 3}, {1, 2, 3}} are not generalized partitions of [3]. 

With the extra restriction that BiDBj = for i ^ j, {Bi, B 2 , ■ ■ ■ , B k } is a partition of 5* 



30| , p. 33]. Recall that the number of such partitions is S(n, k), the Stirling number of the 



second kind, and the number of all partitions equals to B(n), the n-th Bell number. For 
this reason, we also call S(n, k) and B(n) the generalized Stirling numbers of the second 
kind and the generalized Bell numbers, respectively. 

Similar to the recursive formula for Stirling numbers of the second kind, we have the 
following formula for S(n, k). 

Proposition 3.15. For n,k > 1, we have 

(24) 5(n,A;) = 2( n_1 )2^,A;-l). 



S(n,l) = l, S(n,n) = 2&. 

Proof. We first prove the recurrence. 

Every generalized partition of [n] with k blocks is determined uniquely in the following 
three independent steps: 

(1) choose a subset X of [n — 1] with i elements, < i < n — 1; 

(2) choose a generalized partition {£>i, • ■ • , Bk-i} of X; 

(3) choose a subset Bk of [n] that contains ([n — 1]\X) U {n} and contains a subset of 
X. 

The number of choices of the three steps are ("7 1 )' S(i,k — 1) and 2*, respectively, hence 
the recursion. 

Since there is only one way to put [n] in one block, we have S(n, 1) = 1. 

If {B\, • • ■ , B n } is a generalized partition of [n] into n blocks, then we must have max Bi = 
i,l < i < n. Then Bi = {i} U B[ for a subset B[ of [i — 1] (with the convention that [0] = 0. 
Thus there are 2 l ~ l choices for Bi, 1 < % < n. Therefore, we have 

S{n,n) = 2°2 1 ---2 n ~ 1 = 2^). 



It also follows from the recursion. 



□ 
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Proposition 3.16. For n > k > 1, we have 
S(n,k) = 2® 

mi,m 2 ,-" ,m„-k 
1 < mi < m,2 < ■ ■ ■ < m n _k. < n 

Proof. For a generalized partition {B 1: ■ ■ ■ ,B k } of [n] with k blocks. We have max Eli < 
■ ■ ■ < maxBfc. For the convenience of the proof, we use the notation B' maxB . := B i7 1 < % < 
k. In order words, B'- denotes the block whose largest element is j. Such j's form a subset 
N of [n] with k elements. Denote 

[n]\iV = {1 < mi < m>2 < ■ ■ ■ < m n-k < n — 1}- 

With these notations, there are three steps in constructing a generalized partition of [n] 
with k blocks. 

Step 1. Choose a subset iV of [n] of cardinality k to be the largest elements of the blocks. 
Equivalently, choose a subset {1 < m 1 < • • • , m n _fc < n — 1} to be the complement 
of TV. 

Step 2. Determine the contribution of iV in the generalized partition with k blocks. 
Step 3. Determine the contribution of [n]\iV in the generalized partition with k blocks. 

We next derive formulas for Step 2 and Step 3. 
Step 2. We construct a generalized partition of iV with k blocks. We know that all the k 
elements in N are the largest elements of the blocks B' j: j e N. We can also view iV as the 
result after removing the n — k elements of [n]\iV from [n]. Viewing each elements of [n]\iV 
as a cutting point, the elements of N C [n] are cut into n — k + 1 segments even though 
some of the segments might be empty. For example, taking n = 9 and N = {2, 3, 7, 8, 9}. 
Then [9}\N = {mi = 1, m 2 = 4, m 3 = 5, m 4 = 6}, cutting [9] into five groups some of which 
might be empty: 

{}; 

{2,3}; 

{}; 
{}; 

{7,8,9}. 

By collecting the segments together this way, we can put the N elements and their 
corresponding B'- s into groups as follows. 





B' 2 , 


• • * » B ' mi 


TDl 

D mi+H 


TDl 


TDl 


TDl 

- D m 2 + 1' 


TDl 

- D m 2 +2i 


TDl 



TDl TDl . . . TDl 

m„_fc_i+l> - D m, I _ fc _ 1 + l5 i JD m n _ k -l' 

TDl TDl TDl 

- D m n _ fc +1, - D m n _ fc +25 D n- 

We note that B' mi , B' m2 , • • • , B' m _ fe are not defined and don't appear among the above 
groups. Consider Bj with 1 < j < mi — 1, namely the Bj in the first line. Since B'j can 



i=i 

- i 
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contain j and any subset of the j — 1 elements in N that are less than j, there are 2 J_1 choices 
for B'-. Therefore, the total choices of all the blocks in the first line are 2°2 1 ■ • • 2" 11 " 2 . Next 
consider Bj with mi + 1 < j < m 2 — 1. Then Bj contains j and any subset of the elements 
of N that are less than j. There are j — 2 such elements {1, • • • , mi — 1, m\ + 1, ■ • • , j — 1} 
since mi is not in N. Thus there are 2 J_2 choices for Bj. In general, each Bj on the i-th line 
has j — i choices, where 1 < i < n — k+1. Therefore, the number of generalized partitions 
of N with k blocks is. 

~S(N k) = (2°2 1 ■ ■ • 2 mi ~ 2 ) x (2 mi ~ 1 2 mi ■ ■ ■ 2 m2_3 ) x ■ ■ ■ x ^2 mn ~ k ~ n+k 2 mn ~ k ~ n+k+1 ■ ■ ■2 k ~ 1 ) 

= 2 co. 

Step 3: Let a generalized partition of N with k blocks be given. We next determine the 
number of ways to put the elements mi, m 2 , • • ■ , m n _fc into the blocks to form a generalized 
partition of [n] with k blocks. 

By the definition of Bp we see that m.j, 1 < i < n — k can not be put in any Bj with 
j < rrii and hence can only be put in the B'^s after the z-th line in the above table. Since the 
Bj have union [n] and can have overlaps, every m^, 1 < i < n — k can appear in any nonzero 
number of the blocks after the z-th line in the above table. Hence there are 2 k ~( mi ~^ — 1 
possible ways to add to the existing blocks. Therefore, the total number of choices of 
adding [n]\iV to any given generalized partition of iV with k blocks is 

n—k 

^2*-(»ni-i) — 1) x (2 k ~( m2 ~ 2 } — 1) x ■ • ■ x (2 k ~( mn - k ~( n ~ k ^ — 1) = |^J(2 fc ~ m4+i — 1) 

Putting the three steps together, we have 

S(n, k) = 2 ^ n( 2fc " mi+l _ X )' 

mi,m2, ■ ■ • ,m n _ fc 1=1 
1 < mi < r?i2 < ■ ■ ■ < m n -fc < n — I 

We complete the proof. □ 

We have the following special cases of Proposition p. 16 . 
Corollary 3.17. For n > 1, we have 

S{n, 2) = 3™- 1 - 1, S(n, n - 1) = 2^) (2 n -n-l). 

Proof. We give a direct proof of the first equation it can also be obtained from Proposition 
|3.15| and Proposition |3.16| by setting k = 2. 

Let {Bi, B 2 } be a generalized partitions of [n] with 2 blocks. Then 1 < maxi?i < n — 1 
and max B 2 = n. If max Bi = t, 1 < t < n — 1 and j^B\ = j, 1 < j < t, then there are (*Zi) 
choices of B\. For each such a choice of Bi, there are 2 J choices of B 2 since B 2 must contain 
{t + 1, - • • , n} U ([t]\Bi) but can contain any subset of B\. Thus for given 1 < t < n — 1, 
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there are 

E$=i Ql\) 2j = 2(3*- 1 ) choices. Summing over t, we have 

n-l 

S{n,2) = = 3"- 1 - 1. 

i=l 

For S(n,n — 1), setting k = n — 1 in Proposition |3.16| , we have 

n-l 

S{n,n- 1) = 2(V) E (2"- 1 " mi+1 - 1) 



mi=l 
n-l 



2\2 ) {r~ mx - 1) 

mi=l 

2(V)(2 n - n - 1). 



□ 



To obtain the generalization of the Composition Formula, we first need to give a proper 
definition for the composition of two elements in IIIa(IR). We do this by a suitable gener- 
alization of the construction in [25 which treated the case of A = 0. 

Set 



n=0 

For the differential operator cIa, we have ^(ln) = ln-i, and c?a(1) = 0. Then, we have 

oo 
n=0 

Next, for the Rota-Baxter operator Pa, we have P A {l n ) = ln+i- Then we have 

oo 

P A (E XJ ) = J2f(n-l)l n . 

n=l 

Therefore, we have 

d A {P A {E X j)) = E XJ . 

Definition 3.18. We define the n-th divided power E X j[ n ] of E X j recursively by E X j[o\ = 
1 and 

E XJ in) = P A (E XJln -i] o d A (E XJ )) , n > 1. 
Thus for any n > 1, we have g?a(-Ea,/n) = E X j{ n -i\ o dA(E X j). 

Definition 3.19. Given functions / : P — Y R and g : N — > M, we define the composition 

of £ Ai9 and E X j by E x , g (E x j) = Y^=o9i.k)E X f[k] . 

First, we give a formula related to E X j[ k \. 
Lemma 3.20. Given a function f : P — > R, /or k > 1, we have 

oo 

(25) £a,/ W = E E A# Bl+ - + #^-#^ u - u ^)/(#5 1 ) • • • f(#B k )l n , 

n=k {Bi,-,B fc }en'([n]) 
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where {B\, B2 ■ ■ ■ , Bk} ranges over all the generalized partitions of [n]. 
Proof. When k — 1, according to the definition of E X j[k\, we have 

too \ / 00 \ 00 

1 /Hln-1 = P A [J2 f^-l = E / 
n=l / \n=l / n=l 

Setting = 1 in the right hand side of Eq. (|25|) , we have 

00 00 

y: A#*-^/(#soi„=f;/(n)i„. 

™ =1 {Bi}en'([n]) n=i 

Therefore. Eq. ( |~~| ) holds for A; = 1. Assume that Eq. fl25Q holds for fc. Then for k + 1, we 
have 



P AT, E A#^ + - + # B -#^ u - u ^)/(#i? 1 ) • • • 

y k n=fc{Bi,...,B fc }en'([n]) / ^ v v 



o 



\m=l j 



^ EE E A#^+-+# B *-#^ u - u ^)/(#S 1 ) • • • f(#B k )f(m)l n o 1, 

\n=0m=l {Bi,- ,B fc }en'([n]) 

(00 00 
EE E A#^+-+# B -#^ u - u ^)/(#S 1 ) • • • f(#B k )f(m) 

n=0m=l{B 1 ,...,B fc } e n'([n]) 

min(n,m-l) , . . . - 



i=0 

/ 



Pa 



E E 

«=0 ( Ul ,n 2 ,u 3 ) e N 3 

"1 + U2 + «3 = U 



£ A #B 1+ ... + #B fc -#(B 1 U...UB fc)/(#jBi) . . . f{#Bk)f{ui +U3 + 1} ) lw 

{Bi,-,B fe }en'([« 1 +w 2 ]) 



E £ 

u=1 (ux,u 2 ,uz) e N 3 

Ml + U2 + «3 = U — 1 



"1 



u-1 



£ A# B 1+ ... + #B fc -#(B 1 U...UB fc)/(#jBi) . . . mBk)f{ui + U3 + 1} j lu . 

{Bi,-,B fe }en'([«i+« 2 ]) 
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As shown at the end of the proof of Proposition |3.3|, we have 



V" A «i ( ) = V" X #B+#B' k+1 -#(BUB' k+1 ) _ 



(«i,«2,«3)eN 3 x ' {B,B' k+1 }e 7([u-l}) 2 

Ml + U2 + «3 = u — 1 B U B' k+1 = [u — 1] 



Thus we get 



E xjlk+1] =J2 E X #B +m+1 -#(BUB> k+1 ) 



u=1 {B,B' k+1 }en[n-l]) 2 
BUB' k+1 =[u-l] 

• ( A#^ + - + #^-#^ u - u ^)/(#5 1 ) • • • /(#5 fe )/(#^ +1 + 1)) l u . 

{Bi,-,B fe }en'(-B) 

Denote -Bfc+i := U {m}. Then {£>,£>fc + i} forms a generalized partition of the set [u]. 
Further, if {B,Bk+i} is a generalized partition of [u] and {£>i, ■ ■ ■ ,-£>&} is a generalized 
partition of B, then {.B^, ■ • • , -B^+i} is a generalized partition of [u]. Therefore we have 

oo 

E xf[k+1] =^2 E E X #B+#B k+1 -l-(#(BL>B k+1 )-l)+#B 1 + ---+#B k -#B 

u=l {B,B h+1 }£TL'([u}) {B lt - ,B k }€U'(B) 

oo 

= E E A# Bi+ - + #^ + -#^ u - u ^^/(#i? 1 ) • • • 

«=i {Bi,- ,B fc+1 }en'(M) 

oo 

= E E A# Bi+ - + #^-#^ u - u ^)/(#5 1 )---/(#i? fe+1 )i u . 

u =/c+i{B li ... iJ B fc+1 } e n'([«]) 

This completes the induction. □ 

Theorem 3.21. (Generalized Composition Formula) Given functions f : P — > R and 
g : N — > R wrai/i g(0) = I, define a new function h : N — > R fry 

fc(n) = E A^ + - + #^-#^ u - u ^)/(#i? 1 )---/(#i?^(A ; )^ n>0, 

{Bi,-- - ,B fc }en'([n]) 
fc > 1 

h(0) = 1, 

where the sum ranges over all generalized partitions {Bi, • ■ ■ , B k } of [n]. Then 

E\ : h = E\ t g(E\j). 

(Here E\j = J2™=i f( n )ln, since f is only defined on positive integers.) 

Notice that Theorem |3.21| reduces to Theorem [3.13| when A = 0. 
Proof. By Definition [3.19| and Lemma |3.20| , we have 

Ex,g{E X j) 



GENERATING FUNCTIONS AND ROTA-BAXTER ALGEBRAS 



19 



k=0 

oo 

g(0)E XJ[0] +J29(k)E X j [k] 



k=l 

oo oo 



= i+£*(*)£ E A^ + - + ^-#^ u - u ^)/(#fi 1 )---/(#^)in 

k=l n=k {Bi,-,B k }en'([n]) 

oo n 

= 1 + E E E \#^-^-#^-^f^ Bl ) ■ ■ ■ f(#B k )g(k)l n , 

n=l k=l {B u - ,B k }eW([n]) 

as needed. □ 



As applications of Theorem |3.21| , we give the following corollaries, providing further 
justification for the notations S(n, k) and B{n). Recall that E X j = Yl^=o In- 
Corollary 3.22. Let k > 0, E X j = E X j — 1 and E Xi9 = Then we have 

oo 

J2S(n,k)l n , A = 0, 



n=0 

oo 



E X ,g(E X j) = <v 

22S(n,k)l n , A = 1, 

„ n=0 

Especially, when A = 1 and k = 2, we have 



(26) E 1 ,,( J E; 1 , / ) = l + ^(3"- 1 -l)l n . 

71=1 

When A = 1 and k = n — 1, we have 

oo 

(27) E hg (E lif ) = 1 + E 2^ 1 )(2" - n - l)l r 



n=i 



Proof. The first equation follows from Theorem |3.21| since, for our choices, 

/(0) = 0, /(z) = 1, z > 1 and g(i) = 5 kyl , i > 0. 
We then get Eqs. (g§) and (|27|) by Corollary |3T7j . 
By a similar argument, we obtain 



□ 



Corollary 3.23. Let E X j = E X j — 1 — ^2 l n and E Xi9 = E X j = ^2 l n . Then we have 

n=l n=0 

oo 

l + J2 B (n)l n , A = 0, 

n=l 
oo 

l + 5^S(n)l n , A=l, 



E X ,g(E X j) = < 



n=l 
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where B{n) is the n-th Bell number, and B(n) is defined in Definition 14 ■ 



Before proving the next result, we give a definition similar to Definition IXT 



Definition 3.24. For k,£ > 1 and k + £ — 1 < n < k£, let B'(n,k,£) denote the number 
of /c-tuples (Ti, • • • ,Tfc) of (not necessarily disjoint) subsets Ti, • • • , T& of [n] such that 
T\ U • • • U Tfc = [n], #Ti = £ for 1 < i < k, and all the largest elements of the subsets are 
different from each other. For the other values of n, k, £, we set B'(n, k, £) = 0. 

Then by the same proof as for Corollary |3.22| , we have 

Corollary 3.25. If E\j = le an d E\ g = 1^. for k,£>l, then we have 



E Xt9 (E XJ ) = l + J2 B '^ k ^^ 



n=l 

4. Generating functions in free commutative Rota-Baxter algebras 

In this section, we will study generating functions in Ul\(x), the complete free commu- 
tative Rota-Baxter algebra on one generator x. In Section [ll] and jOj we will discuss the 
weight zero case and the general weight case, respectively. We first recall the notations: 

00 

HI A (x) := ni A (M[x]) = R[xf k = R{x no <g> x ni <g> ■ ■ ■ ® x nk \ n Q , ■ ■ ■ , n k > 0, k > 0}. 

k=l 

Recall that the product in HIa(x) and hence in IU\(x) is given by the mixable shuffle 
product o which will be suppressed. 

We state some definitions and notations which are used in this section. A composition 
of n is an expression of n as an ordered sum of positive integers. If exactly k summands 
appear in a composition a, we say that a has k parts, and we call a a ^-composition. Let 
/ = (ii, ■ ■ ■ ,i t ) G P*. Define the norm of / to be |/| = J^* =1 i s and its length to be 

n \ / n 

I I \ ii,— ,it 



£{I) = t. Also denote 



4.1. The weight zero case. We consider the case when A = 0. Then for x n ° (g) • • • <g) x Uk 

and x m ° (g) - • • (g) x mt we have 

(X n ° ® • •• ®X nk )(x m ° ®---®X mt ) = X n0+m ° ® ((X ni (g) • • • g> X^m^x" 11 g) •• • (g X mi )), 

where m is the shuffle product. In particular, 

(1 <g> x® n )(l ® x® m ) = ( U t m ) 1 ® x^ (n+m) 

and more generally, 



(28) TT(l<g^) = ( Ul + '" + nk )(l®x ri 

Also, 



(niH hn k )\ 



(29) n( l0a;i ) = Yl 1 ® xa(1) ■ 



i=i o-es k 
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where Sk is the set of permutations on [A;]. 
Then we easily obtain 

Proposition 4.1. For n,k>l,we have the following relations in ILT (a;). 
(30) 



ra / (n+l\ \ 

rr ( i® x ^)= u J (i®^ 1 )). 

fc=1 V 1 ; A 3, . . . , n/ 



Further, 



n=l n=0 |J|=n 

where J denotes all the compositions of n with distinct parts. 



Proof. The first two equations follow from Eq. (|28|) . The third equation follows from 
Eq. (HD. □ 

By treating 1 ®x® k G \Hq{X) as the base q in basic hypergeometric series, we next apply 
g-series identities to derive generating functions for some special sequences. 



Recall the g-shifted factorial \. 



' ?)°° = TT(1 - aq n and (a; q) n = °* 9 °° , for n € Z. 
1 ' (aq n ;q) OQ 



(a; g) 

l^i - ag ; anu i^u; g; n = - 
n=l ' ' 

Then we have the following identities || in 

oo 

(31) £ g n2 = (-g;g 2 ) 2 oo(? 2 ;? 2 )oo, 

n=— oo 



oo 



(32) *h):=E ' ^ 



n=0 

oo 



n(3n — 1) 



2^ 

oo 



(33) /(-g) := £ ("1)V^ = (9; 9)oo. 

n=— 00 

We have the following generating functions for the sequences of the factorial of square 
numbers, triangular numbers, and pentagonal numbers. 

Theorem 4.2. We have 

00 00 

(™ 2 )'(1 ® ^™ 2 ) = II^ 1 + ( 2n ~ W ® ^ (2n_1) )) 2 ( 1 - (2n)!(l ® a;® 2n )), 

n=— 00 n=l 

TT (l-(2n)!(l®x^ 2n ))) 



^\ 2 j' 1 ^ l\ (l-(2n- 1)!(1 ®^(*-D))' 

E ( n(3n 2 " 1) )l(l ® s* 2 ^) = f[(l - n!(l ® *«»)). 
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Proof. Note that M [[(?]] is the free completed commutative K-algebra on q in the sense that 
for any complete IR-algebra A and nonunit a G A, there is unique algebra homomorphism 
g : M.[[q\] — > A such that g(q) = a. Thus about identities in M.[[q]] induce similar identities 
in A. Take A = IIl(x), a = l®x and let g x : M[[g]] — > Ul(x) be the algebra homomorphism 
such that g x (q) — 1 <8> x. Applying g x to the above three g-series identities in M[[g]], we 
obtain identities for x G m(x). 

For example, applying g x to Eq. (|31~D, we have 

oo / oo \ ^ / oo 

J2 (1 8) a;)" 2 = m (1 + (1 <g) a;f t 2 "- 1 )) I JJ(1 - (1 <g> x) 2n " 

n=— oo \n=l / \n=l 

Then using the following special case of (|30| ) 

(1 ®x) n = n\(l®x® n ), 
we obtain the first equation. The proofs of the other two equations are similar. □ 

4.2. The general weight case. Before giving the main theorem in this section, we first 
need some definitions and notations. 

Let ir(kn) denote the set of compositions of kn with each part less or equal to n. For a 
composition / = ^2? " £ ^{kn), we define a restricted partition of type / of 

the multiset S := S n ^ '■= {l k , 2 fc , • • • , n k } to be a partition {Bi, • • • , £? t } of 5 such that 

(1) i,. 1 < ./ < /: 

(2) each 5j is a subset of [n). In other words, #(Bj fl {j fe }) < 1; 

Let Cj := C nt kj denote the number of the restricted partitions of S of type /. We will 
use the convention that, for 7 = 0, take I G 7r(0), Cj = 1, £(I) = kn, and 1 (g> = 1. 



Proposition 4.3. For n, /c > 1, we /iave 

(kn)\ 
= (k\) n 





and 



Proof. For I G 7r(A;n), if I = ( 1,1,-- ,1 ), and the multiset S = {l k , 2 k , ■ ■ ■ , n k }, Then Cj 

denotes the number of ways to put the kn elements in S into kn different subsets, where 
each subset contains 1 element. Therefore, we have 

(kn)\ 
= (k\) n ' 

If I — (n,n, ■ • ■ , n) and S = {l k , 2 k , ■ ■ ■ , n k }, then Cj denotes the number of ways to 





put the kn elements in S into k different subsets, where each subset contains n different 
elements. There is only one way to do this, hence the second equation. □ 

Proposition 4.4. For n, k > 1, we have 

Cj = C(n } k) } 

/G7r(fcn) 



GENERATING FUNCTIONS AND ROTA-BAXTER ALGEBRAS 23 

where C(n,k) denotes the number of ways to partition kn elements in a multiset S = 
{l k , 2 k , ■ ■ ■ , n k } into several nonempty subsets of [n]. 

Proof. For all the compositions in 7c(kn), according to the definition of Cj, we don't need 
to limit the number of elements in each subset of [n]. □ 

Tracing back to Cartier's work [7j forty years ago, we also recall the description of the 
mixable shuffle product in terms of order preserving maps ]17], § 3.1.4] that are called 
stuffles in the study of multiple zeta values |5|, H, HD- We will use this interpretation in 
the proof of Theorem [L5j 

For < r < rninfm, n), define 



Jm,n,r = \ V>) 



ip : [m] — > [m + n — r] , ip : [n] — > [m + n — r] are order preserving 
injective maps and im(ip) U im(ip) = [m + n — r] 

Let a G A® m , b G A® n and (</?, ip) G J m ,n,r- Define ami^^b to be the pure tensor of length 
m + n — r whose i-th factor is 

{dj, if i = tp(j), i G" fnT0, 
bj, if i = ip(j), i £ imy?, 
ajb jl} if i = <p(j), i = ip(f). 

In short, 

a^{ip,4>)b = CL v -i{i)b^-i{i) ® ■ ■ ■ <8> a^-i( TO+n _ r )&^-i( m+n _ r ), 
with the convention that a% = b% = 1. Then dm(„^6 is called a stuffie of a and 6. The 
stuffle product of a and b is defined by 

min(m,n) / 

am st,\b — A r I Om(^)6 

r=0 \(p,lt)6Jm,n,r 



which coincides with the mixable shuffle product 
Theorem 4.5. We have the generating function 



1 - 1 ® x 0/£ 

n=0 n=0 l£ir{kn) 

Proof. According to the description of the mixable shuffle product in terms of stuffles stated 
at the beginning of this section, for (1 £g> x® k ) n , we define an n-fold stuffle to be 

fi : [k] — > [t] for 1 < % < n are order preserving 
injective maps and U™ =1 im fi = [t] 



£>k,n,t — { {/lj /2; " " " , fn} 

Set 



kn 



£>k,n,t — £>k,n,t- 



t=h 

„®k\n 



Then for each {/i, / 2 , ■ ■ • , /„} G £>k,n,t, we get a term in the expansion of (1 £g> x® ) n as 

1 <g> x' 1 <g> x i2 <g> • • ■ ® a:**, 



where 



i = #(U£=i/- 1 (j))forl<i<*- 
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Further the following relations hold. 



t 



1 < ij < n for 1 < j < t and ij = kn. 

Therefore, each term of (1 <S> x® k ) n has the form 1 eg) x® 1 , where / is a composition of kn 
with each part less or equal to n. 

Given I — (ii, ■ ■ • , i t ) e 7r(kn), we also define an n-io\d stuffie of type / to be 

fi : [k] — >■ [i] for 1 < i < n are order preserving injective maps, 
such that U™ =1 im/< = [t] and #(U? =1 / _1 (j)) = ij, l<j<t. 

Let 5"/ denote the number of n-fold stuffles of type /. Then we have 



{fl, f2, ■ ■ ■ , fn} 



n=0 n=0 l£ir(kn) 

Thus we just need to prove Cj = Sj. 

Let {fi : [k] — > [t]}i<i< n be a stuflie of type /. We can derive from /j a 

g:[n)x [k] -> [t], g(i,j) = g^j), G [n] x [k). 

Identifying [n] x [k] with the multiset S = {l k , 2 k , ■ ■ ■ , n k }, we obtain a restricted partition 
of 5* of type I by defining 

B r .= g- l {j),l<j<t. 

Conversely, given a restricted partition {-Bj}i<j<t of 5 of type I and identifying it with 
an ordered partition of [n] x [k]. Then define 

g:=g B : [n] x [fc] ^ [t], j) = £ if (i, j) e B t . 

Then it is easy to check that the maps 

fi : [k] -)> [t], := ff(z,j), 1 < 2 < n, 

is a stuflie of type /. This gives the desired bijection that gives Ci = Si. □ 



Based on the definition of Cj and the proof of Theorem |4.5| , we have the following result. 
Corollary 4.6. Forn,k > 1, we have 



kn 



J2B(t,n,k) = C(n } k). 



t=k 



Proof. Given ale ir(kn), in the proof of Theorem [43], we know that Ci = Si, where Si 
denotes the number of n-fold stuffles of type /. Therefore, we have a equivalent illustration 
for Ci. For any composition I = (i 1 ,i 2 , ■ ■ ■ , U) G iz{kn), Ci denotes the number of n-tuples 
(Ti, T 2 , ■ • • , T n ) of subsets T\, T 2 , • • • ,T n of [t] such that = k and all the elements in 
Ti,T 2 , ••• ,T n are composed by ij j's, for 1 < j < t. If we take all the t-compositions 
in ix(kn), then we don't need to consider the number of appearances of j for 1 < j < t. 
Therefore, we have 

Ci = B{t,n,k), 

I g n(kn) 
<(/) = t 
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where B(t,n, k) is explained in Definition p.llj . According to Proposition |4.4| , we complete 
the proof. □ 



Now we consider some special cases of Theorem [4.5| for A = 0, 1. For A = 0, we have the 
following generating function. 

Corollary 4.7. We have 

- £(1 ® = £ ® G mo(x). 



l-l<g>x® fc ^ ^-Wfc!) n 

n=0 n=0 V ' 



Proof. When A = 0, Theorem [15] reduces to 

1 oo oo 

!_ 18IJ » =E £ C,(l«*«) = 5: 1,(1 (8 ****>) 

«=0 7 g Tr(fen) "=0 S ~j£ ' 

<!(/) = kn 

which gives what we need by Proposition [O. □ 



Setting x = 1 in Corollary |4.7| , we obtain Eq. 
By taking the special case when A = 1 and k — 1 in Theorem [4.5| , we obtain the following 
result in [16[ . 

Corollary 4.8. We have the generating function 



1 



,1 « x) 



OO OO /• \ 

n=0 n=0 in=n ^ ' 



where for 7 = 0, ia£;e |/| = 0, ^'J'j = 1, (1 <g> x) J = 1 g> x® 7 = 1. 
Proof. Setting A = 1 and k — 1 in Theorem [O, we have 



1 - 1 (8> x 

n=0 n=0 Jg7r(n) 

where ?r(n) is the set of all the compositions of n. For a composition J = (ii, ■ ■ ■ , i t ) G vr(n), 
the coefficient Cj is the number of ways to divide the n elements in [n] into k disjoint 
nonempty subsets, where the j-th subset contains ij elements, 1 < j < k. So we get 

Ct = ( n I , as needed. □ 
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